A Kinematically Complete Analysis of the CLAS Data on the Proton Structure 

Function F 2 in a Regge-Dual Model 



R. Fiore Q , A. Flachi fc , L.L. Jenkovszky c , A.I. Lengyel d , and V.K. Magas c ' e 
a Dipartimento di Fisica, Universita della Calabria & Instituto Nazionale di Fisica Nucleare, 
Gruppo collegato di Cosenza, 1-87036 Arcavacata di Rende, Cosenza, Italy 
b IFAE, Campus UAB, 08193 Bellaterra (Barcelona), Spain 
c Bogolyubov Institute for Theoretical Physics, Academy of Sciences of Ukraine, UA-03143 Kiev, Ukraine 
d Institute of Electron Physics, Universitetska 21, UA-88000 Uzhgorod, Ukraine 
e Center for Physics of Fundamental Interactions ( CFIF), and 
Physics Department, Instituto Superior Tecnico, Av. Rovisco Pais, 1049-001 Lisbon, Portugal 

(Dated: February 1, 2008) 

The recently measured inclusive electron-proton cross section in the nucleon resonance region, 
performed with the CLAS detector at the Thomas Jefferson Laboratory, has provided new data 
for the nucleon structure function F2 with previously unavailable precision. In this paper we 
propose a description of these experimental data based on a Regge-dual model for F2. The 
basic inputs in the model are nonlinear complex Regge trajectories producing both isobar 
resonances and a smooth background. The model is tested against the experimental data, and 
the Q 2 — dependence of the moments is calculated. The fitted model for the structure function 
(inclusive cross section) is a limiting case of the more general scattering amplitude equally ap- 
plicable to deeply virtual Compton scattering (DVCS). The connection between the two is discussed. 
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PACS numbers: 12.40.Nn, 13.60.Hb 14.20.Dh 
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00 . I. INTRODUCTION 

i> 

qq ' It has been recently realized QQHiBI that a straightforward generalization of the ordinary parton densities arises 
in exclusive two-photon processes in the so-called generalized Bjorken region, e.g. in Compton scattering with a highly 
C*~) 1 virtual incoming photon, and in the hard photoproduction of mesons. Here one finds off-forward matrix elements, as 
I distinguished from the forward ones in inclusive reactions. 

Deeply virtual Compton scattering (DVCS) combines the features of the inelastic processes with those of an elastic 
Q_i| process. The diagram of such a process, e(ki) +p\ — ► e'{k 2 ) + P2 + 7(92), is shown in Fig. ^ where e(fci), e'(k 2 ) 
JL . denote, respectively, the initial and final electrons of momenta k\, k 2l andpi, p 2 denote the initial and final momenta 
of the target correspondingly. 

> 
• 1— 1 

x 




FIG. 1: Kinematic of deeply virtual Compton scattering. 

DVCS is the hard electro-production of a real photon, i.e. 7*iV — * jN', Being a process involving a single hadron, 
it is one of the cleanest tools to construct generalized parton distributions (GPD) [E El IS @ i which reduce to 
ordinary parton distributions in the forward direction. The theoretical efforts and achievements are supported by the 
experimental results from HERMES, HERA and CLAS Collaborations, and encouraging future plans. 
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DVCS is characterized by three independent four-momenta: p = p\ + P2, A = pi — p\, and q = (qi + ?2)/2, where 
the vectors pi (qx) and P2 (92) refer to the incoming and outgoing proton (photon) momentum, respectively. Most of 
the papers on deep inelastic scattering (DIS) and DVCS are based on the operator product expansion with extensive 
use of the light-front variables. Otherwise, the conventional Bjorken variable is x = 2 p 1 . qi 1 Q 2 = ~Qi> an d £ = — 
is the generalized Bjorken variable. If both photons were virtual, we would have an extra scaling variable r\ = the 
skewedness (or skewness) |TL lldj) . The reality of the outgoing photon implies the presence of only one scaling variable, 
namely, for q\ = one has 



A 2 \-i 
2Q 2 



V = -^ + 7^) ■ (1) 



The generalized and ordinary Bjorken variables are related by 

A 

f 

a; A, 



£ = * 1 ' 2 °V ( 2 ) 



Our starting point is a complex scattering amplitude depending on three variables, £,i and Q 2 , defined by Fig. 1 
and the corresponding legend. Even though our paper is devoted to DIS of Fig. 2 and relevant CLAS data, we bear 
in mind the close relation between DVCS and DIS, the latter being the limiting case of the former. 

Most of the papers on this subject are based on the factorization properties, separating the perturbative and non- 
perturbative dynamics ("handbag" diagram), according to which, at large Q 2 , lowest order perturbation decouples 
from hadronic dynamics during the short time of interaction. While factorization in hard scattering processes is valid 
to all orders in perturbation theory, a considerable fraction of the existing data comes from so-called soft region of 
small and intermediate values of Q 2 (Q 2 ~ 1 GeV 2 ), where the present non-perturbative approach can be compared 
with the relevant successful pQCD calculations 0, [| ■ Although t— dependence at small t is outside the perturbative 
QCD (pQCD) domain, nontrivial forms of the t— dependence at a proper scale suggested recently 0, 0, 0] can be 
confronted to those following from Regge-dual models. 

The phase of the DVCS amplitude experimentally is extracted from the interference between the DVCS and Bethe- 
Heitler amplitudes, like in the case of the Coulomb interference in the forward cone of elastic hadron scattering. While 
pQCD factorization details 0, 13 how to calculate the real part of the DVCS amplitude, any Regge-dual model 
contains the phase explicitly, its form depending on the available freedom (form of the Regge singularity, shape of the 
trajectories etc) inherent in this type of models. One can hope that the results of the pQCD calculation will reduce 
this freedom in the future. Alternatively, this phase can be approximately reconstructed by means of the dispersion 
relations or their simplified version of the derivative dispersion relations, as it was done in ref. 

In a series of papers we initiated the study of DIS and DVCS within a Regge-dual approach. Its virtue is the 
presence in the scattering amplitude of t— dependence and of the phase as well as its explicit energy dependence, 
compatible with unitarity. At high energies, the contribution of a dipole pomeron [T^ dominates, while at moderate 
and low energies subleading contributions (secondary reggeons) become important. Moreover, by duality, at low 
energies, t— channel Regge pole exchanges are replaced by direct-channel reggeons. 

No hard scale factorization is assumed in this approach. External photons interact with the proton via vector meson 
(or generalized vector meson 

) 

dominance. 

The main idea behind the model is reggeization of the resonances both in the s— and t— channels. Nonlinear, 
complex Regge trajectories replace individual resonance contributions. The resulting scattering amplitude is a complex 
function of the Mandelstam variable s,t,u and of the photon virtuality Q 2 . Its imaginary part in the forward 
direction, t — corresponds to ordinary distributions or structure functions (SF), describing inclusive (e.g. electron- 
proton) scattering, while the whole amplitude is directly related to exclusive deeply virtual Compton scattering and 
corresponding gen eral parton distributions. 

In Refs. [Ij, [IfJ [13 dual amplitudes with Mandelstam analyticity (DAMA) were suggested as a model for DVCS 
or DIS. We remind that DAMA realizes duality between direct-channel resonances and high-energy Regge behavior 
( " Veneziano-duality" ) . By introducing Q 2 -dependence in DAMA, we have extended the model off mass shell and 
have shown 0, now parton-hadron (or "Bloom-Gilman") duality is realized in this model. With the above 
specification, DAMA can serve as explicit model valid, in principle, at all values of the Mandelstam variables s, t and 
u as well as for any Q 2 , thus realizing duality "in two dimensions": between hadrons and partons, on the one hand 
and between resonances and Regge behavior, on the other hand. The latter property opens the way of linking JLab 
(large x, resonances) and HERA (small x, Regge) physics. 

Recently new data on inclusive electron-proton cross section in the resonance region (W < 2.5 GeV) at momentum 
transfers Q 2 below 4.7 (GeV/c) 2 , measured at the JLab (CEBAF) with the CLAS detector 01 were made public. In 
the present paper we discuss an analysis of the new CLAS data within this model. 
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The kinematics of inclusive electron-nucleon scattering, applicable to both high energies, typical of HERA, and low 
energies as at JLab, is shown in Fig. [2] (see ref. ^(| for more details). 




FIG. 2: Kinematic of deep inelastic scattering. 

Studies of the complex pattern of the nucleon structure function in the resonance region have a long history (see, 
for example Among dozens of resonances in the j*p system above the pion-nucleon threshold only a few of 

them can be identified more or less unambiguously. Therefore, instead of identifying each resonance, one considers 
a few maxima above the elastic scattering peak, corresponding to some "effective" resonance contributions. Recent 
results from the JLab 0,^1 renewed the interest in the subject and they call for a more detailed phenomenological 
analysis of the data and a better understanding of the underlying dynamics. 

The basic idea in our approach is the use the off-mass-shell continuation of the dual amplitude with nonlinear 
complex Regge trajectories. We adopt the two-component picture of strong interactions, according to which direct- 
channel resonances are dual to cross-channel Regge exchanges and the smooth background in the s— channel is dual 
to the Pomeron exchange in the t— channel. As explained in ref. |14| . the background in dual model corresponds to 
pole terms with exotic trajectories that do not produce any resonance. 



II. REGGE-DUAL STRUCTURE FUNCTION 



In the present section we introduce notations, kinematics and the Regge-dual model. More details on the model 
can be found in earlier paper [l4l ITEl Il6l [2(j ■ 

So, we study inclusive, inelastic electron-proton scattering, whose cross section was measured at JLab and used to 
determine the unpolarized structure function £2(0;, Q 2 ) as well as the Nachtmann and Cornwall-Norton moments (see 
e.g. 0). 

The cross section is related to the structure function by 

where the total cross section, a] p , includes by unitarity all possible intermediate states allowed by energy and 
quantum number conservation, and we follow the norm 

af p (s)=lm A(s,Q 2 ) . (4) 

used in Refs. 0, flit liH l22j . The center of mass energy of the j*p system, the negative squared photon virtuality 
Q 2 and the Bjorken variable x are related by 

s = W 2 = Q 2 (l-x)/x + m 2 . (5) 

In the Regge-dual approach with vector meson dominance implied, Compton scattering can be viewed as an off- 
mass shell continuation of a hadronic reaction, dominated in the resonance region by non-strange (N and A) baryonic 
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resonances. The scattering amplitude can be written as a pole decomposition of the dual amplitude and factorizes as 
a product of two vertices (form factors) times the propagator: 

where N is an overall normalization coefficient, r runs over all trajectories allowed by quantum number conservation 
(in our case r = N*, N$, A) while n runs from n™ m (spin of the first resonance) to n™ ax (spin of the last resonance 
- for more details see next section), and [A(s, Q 2 )]^o is the contribution from the background. The functions f r (Q 2 ) 
and a r (s) are respectively form factors and Regge trajectory corresponding to the r*' 1 — term. (For a comparison of 
the direct-channel, "reggeized" formula J^J with the usual Breit-Wigner expression see Appendix 0. Note that only 
for the first resonance at each trajectory we have squared form factor, while for the recurrences the powers of form 
factors are growing, according to the properties of DAMA 0, . 



A. Regge Trajectories 



Any systematic account for the large number of direct-channel resonances (over 20) contributing to the 7* N total 
cross section with different weights is not an easy task. However, this problem can be overcome with the use of (s- 
channel) Regge trajectories, including all possible intermediate states in the resonance region appearing as recurrences 
on the trajectories. In this approach, Regge trajectories play the role of dynamical variables and the parameters of 
the trajectories can be fitted either to the masses and widths of the known resonances or to the data on DIS cross 
sections (structure functions), reflecting adequately the position of the peaks in the SF (or cross sections) formed by 
the interplay of different resonances. 

The form of the Regge trajectories is constrained by analyticity, requiring the presence of threshold singularities, 
and by their asymptotic behaviour imposing an upper bound on their real part. Explicit models of Regge trajectories 
realizing these requirement were studied in a number of papers |23| . For our present goals (small- and intermediate 
energies) a particularly simple model based on a sum of square root thresholds will be suitable: 

a(s) = a Q + ats + a 2 (V^o _ V«o - s), (7) 

where the lightest threshold, sq, produces the imaginary part and the heaves thresholds producing the real part can 
be approximated here by a linear term. In our case |l4Lll5lll6|| sq = (m v + m p ) 2 . 

For asymptotic, large s— or t— the trajectories turn down to a logarithm, producing wide angle scaling behavior 
with a link to the quark model. This interesting regime, discussed e.g. in ref. however is far away from the 

resonance region and will not be included in the present analyses. 

In 7*p scattering, mainly the two A*s (isospin 1/2) and one A (isospin 3/2) resonances contribute in the s— channel 
and thus we will limit ourselves to considering these three terms, plus an additional terms which describe the back- 
ground, to be be discussed later. 



B. Form Factors 

In our previous work [l6|. we concentrated our attention on the analytic structure of the scattering amplitudes 
using a simple dipolc model for the form factors. However, in order to properly describe the structure function in 
the resonance region, it is essential to account for the helicity structure of the amplitudes. Below we do so following 
Davidovsky and Struminsky [26j, who provided for relevant amplitudes by using the Breit-Wigner resonance model. 
The relation between the Breit-Wigner and the "reggeized" resonance model, to be used can be found in Appendix 

El 

The form factors can be written as a sum of three terms j2|| I27I I2H . Efll ] , G+ (Q 2 ), Go(Q 2 ) and G_ (Q 2 ) , corresponding 
to 7* A — > R helicity transition amplitudes in the rest frame of the resonance R: 

<R,X R =X N -X 1 \J(0)\N,X N > 
m 

where Xr, Xn and A 7 are the resonance, nucleon and photon helicities, J(0) is the current operator; A 7 takes the 
values —1, and +1. Correspondingly, the squared form factor is given by a sum |2fil 1271 l28t l29j 

|G+(Q 2 )| 2 + 2|Go(Q 2 )| 2 + |G_(Q 2 )| 2 . (9) 
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The explicit form of these form factors is known only near their thresholds \q\ — > 0, while their large-Q 2 behavior 
is constrained b y th e quark counting rules. 
According to '2T\ , one has near the threshold 



G ± {Q 2 )~\q\ J -"\ G (Q 2 ) ~ |kY- 1/2 
for the so-called normal (l/2 + — * 3/2~, 5/2 + , 7/2~, ...) transitions and 

G±(Q 2 ~ W- 1 ' 2 , G Q (Q 2 ) ~ |k1 J+1/2 
for the anomalous (l/2 + — > l/2~, 3/2 + ,5/2~,...) transitions, where 



(10) 



(11) 



kil 



yj (M 2 - m 2 - Q 2 ) 2 + 4M 2 Q 2 
2M 



(12) 



M 2 - m 2 



9o = 



2M 



(13) 



M is a resonance mass. 

Following the quark counting rules, in refs. [2!| (for a recent treatment see |2(|), the large-Q 2 behavior of G's was 
assumed to be 



G+ ~ Q 3 , Gq ~ Q 4 , G_ ~ Q 



(14) 



Let us note that while this is reasonable (modulo logarithmic factors) for elastic form factors, it may not be true any 
more for inelastic (transition) form factors. Our Regge-dual model, eq. JJjJ, predicts that the powers of the form 
factors increase with increasing excitation (resonance spin). This discrepancy can be resolved only experimentally, 
although a model-independent analysis of the Q 2 -dependence for various nuclear excitations is biased by the (unknown) 
background. 

In ref. [2(| the following expressions for the G's, combining the above threshold- I|1U|) . I|ll|l with the asymptotic 
behavior Q14p. were suggested: 



|G±| 2 = |G±(0)| 2 q 



2 „2J-3 „2J-3( n /\ _m± 



(Go) c m± (Q ) 



(15) 



for the normal transitions and 



|Go| 2 = G 2 ^g 2J - 1 C 2a+ " l «(Qo)c 2J - 1 (Q ) 



(16) 



|G±| 2 = \G ± (0)\ 2 q 2J - 1 c 2J - l (Q' ) c m ±(Q ) 



(17) 



l G o| 2 = c2 (|^) 2J 1 c 2a+m °(Qo)c 2J+1 (Q' ) , (18) 

for the anomalous ones, where m + — 3, mo = 4, m_ = 5 count the quarks, G and a are free parameters. For 
notational convenience we have introduced the functions 



\o\q=o ' 

,2 



) 2 + Z 2 



The form factors at Q = are related to the helicity photoproduction amplitudes Ai/ 2 and A 3 / 2 by 



|G+,_(0)| - -*J-¥—\A 1/2 , a/2 \ . (19) 
V Ana \ M —m 
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C. The Background 



Apart from the resonances, lying on the iV"*'s and A s— channel trajectories, dual to an effective bosonic 
(/-) trajectory in the t— channel, one has to consider the contribution from a smooth background. Following our 
previous arguments, 

CI 13 El IH, we model it by non-resonance pole terms with exotic trajectories, dual to the 

Pomeron. 

[M^Q 2 )] BG = £ G » (2°) 

b=E,E' n b- a b{ s ) 

with dipole form factors, given by c 2 (Qb)- The exotic trajectories are chosen in the form 

a b (s) = a<6(0) + aibiy^ - \/s - s) , (21) 

where the coefficients a&(0), a\b and the Q 2 are the free parameters. To prevent any physical resonance, they are 
constrained in such a way that the the real part of the trajectory terminates before reaching the first resonance on the 
physical sheet. An infinite sequence of poles, saturating duality, appears on the non-physical sheet in the amplitude; 
they do not interfere in the smooth behaviour of the background (for more details see [20]). 

Anticipating the results of Sec. IV, we notice that fits to the data prefer a ne gat ive contribution from the second 
term in the background. Formally this is compatible with alternative models (e.g. |17L fl9| 't ■ but needs to be understood 
also in the framework of the present Regge-dual approach. 



III. COMPARISON WITH OTHER MODELS 



In this section we would like to indicate on the two important properties of our Regge-dual model, that should, in 
principle, discriminate it from alternative models of DIS in the resonance region. 

Looking at eq. JfjJ) one can see that contrary to the models accounting for each resonance separately here resonances 
on each Regge trajectory enter with progressively increasing powers of the form factors. This makes the present model 
quite different from the existing approaches [l^.ll8LIT9ll2r| . Notice that increasing powers of the transition form factors 
result in the suppression of the relevant contributions from the recurrences with growing spin, thus explaining the 
gradual disappearance of higher excitations. Further comparison of the experimentally measured transition form 
factors may discriminate between two approaches. Work in this direction is in progress. 

The second important difference comes form the parameterization of the background. We describe the background 
by non-resonating pole terms (the poles appear on the non-physical sheet, see |30j) with exotic trajectories and stan- 
dard dipole form factors. The background contribution strongly decreases with increasing Q 2 , whereas in "standard" 
parameterizations 0,0,^1113 the background is an increasing function of Q 2 . Since resulting fits by different mod- 
els are almost equally good, it is difficult to discriminate between these two options. Studies of the Q 2 — dependence 
of the ratio between a resonance contributions and the background (at fix energy or x) may resolve this ambiguity 
and help to better disentangle resonances from the background. 



IV. ANALYSIS OF THE CLAS DATA 

In this section, we present our fits to the CLAS data on the nucleon structure function, F^Q 2 , s) 
A similar analysis using earlier data |19| was carried out in our previous paper [lq|. The main point of the 
model considered in [16| was the inclusion of three prominent resonances, iV*(1520), A r *(1680) and A(1232) plus a 
background, dual to the Pomeron exchange. In that approach the large number of resonances contributing to the 
i*2 with different wights was effectively accounted for by letting the SF to depend on effective trajectories, whose 
parameters were fitted to the data. This approach was, in a sense, justified "a posteriori" : the parameters of the 
effective trajectories were found to be close to these fitted to the spectrum of baryon resonance. Although the main 
features of the SF in 0] were reproduced by the dual model, the quality of the fit was far from perfect. The reason 
for the poor agreement could be threefold: first, in EU we made an extra simplification by neglecting the helicity 
structure of the amplitudes, and the form factors were chosen in a simple dipole form. Including the spin changes 
the form factors in a non-trivial way and complicates the Q 2 — dependence of the SF. The second point is related 
to the parameterization of the background: in [l(| the background was modeled by one term only, underestimating 
the magnitude of the SF in some regions. The third important reason is the quality of the data - the set of points 
available was not homogeneous resulting in a non- uniform weight of the fit. To cure this deficiency, we performed a 
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preselection of the initial data set, a procedure that potentially may results in ambiguities. The fits were improved, 
although still are not perfect. 
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FIG. 3: Comparison between three different fits performed in the present model (see text). 



Similarly to [16(, here we also include only the contribution from three dominant resonances: AT* (1520), TV* (1680) 
and A(1232) and we implement this by using three baryon trajectories with one resonance on each of them. By 
considering such resonances as "effective" contributions to the SF, we are able to treat the large number of resonances 
that contribute, with different weights, to the SF. 

The imaginary part of the scattering amplitude can then be written, according to © , as a sum of the contribution 
from the resonances plus the background, 

Im A(s, Q 2 ) = N{ [Im A(s, Q 2 )] R + [Im A(s, Q 2 )]bg) ■ 

Accordingly, the resonance contribution takes the following form: 



[Im A(s,Q 2 )] R = Yl fHQ 2 

j=A, N u N 2 



Im, 



(rij - TZej) 2 +lm 2 j ' 



with IZe and Im denoting the real and imaginary part of the relevant Regge trajectory, and the form factors are 
calculated as described in sec. Ill Bl For instance, the form factor for the A resonance can be written as 



fl(Q 2 ) = q 2 c 2 (Q' ) (c 3 (Q )|G + (0)| 2 + c 5 (Q )|G_(0)| 2 



(22) 
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similar expressions can be cast for other contributions. 

The imaginary part of the forward scattering amplitude coming from the background can be easily obtained from 



[1m A{ Sl Q 2 )] BG 



E 

j=E,E> 



Im, 



(nf in -Kej) 2 +lm 2 



Here n™ m is the lowest integer, larger than Max [Hej] , ensuring that no resonances will appear on the exotic trajectory. 
The advantage of such choice is that the two terms of the background depend on two different scales, Q% and Q 2 E , , 
so they will dominate in different regions. 

The model constructed in this way, has 23 free parameters: each resonance is characterized by three (the intercept 
is kept fixed) coefficients describing the relevant Regge trajectory plus the two helicity photoproduction amplitudes 
(see eq. O)- The form factors (see sec. Ill Bfl leave only two free parameters, Qo an d Q'g- Finally, the background, 
contains 8 free parameters: 4 for the two exotic trajectories, 2 energy scales Qe and Qe> and two amplitudes Ge and 
Ge 1 - With the overall normalization factor, N this gives a total of 23 free parameters. 

The resulting fits to the CLAS data, performed by using MINUIT [3l|. are presented in Table 1 and together with 
the experimental data is shown for various Q 2 bins in Figs. 1)71 12[l . 
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ao 




-0.8377° 


-0.8377° 


-0.8377° 




ai [GeV" 2 ] 




0.9500° 


0.9402 


0.9825 




a 2 [GeV" 1 ] 




0.1473° 


0.1757 


0.0920 




A 2 (1/2) [GeV" 




0.0484E-2° 


0.0484E-2° 


0.864 7E-2 




A 2 (3/2) [GeV~ 


"] 


0.2789E-1° 


0.2789E-1° 


0.9634E-2 




QO 




-0.3700° 


-0.3700° 


-0.3700° 




ai [GeV" 2 ] 




0.9500° 


0.9724 


0.9551 


iV| 


a 2 [GeV" 1 ] 




0.1471° 


0.0575 


0.0949 




A 2 (1/2) [GeV" 


'} 


0.0289E-2° 


0.0289E-2° 


0.9724E-2 




A 2 (3/2) [GeV" 


'} 


0.1613° 


0.1613° 


5.1973E-11 




QO 




0.0038° 


0.0038° 


0.0038° 




ai [GeV" 2 ] 




0.8500° 


0.8758 


0.8605 


A 


a 2 [GeV" 1 ] 




0.1969° 


0.1724 


0.2005 




A 2 (1/2) [GeV" 


') 


0.0199° 


0.0199° 


5.3432E-08 




A 2 (3/2) [GeV" 


'] 


0.0666° 


0.0666° 


0.0866 




Ge 1 




6.5488 


2.8473 


3.6049 




ao 




0.3635 


0.7014 


0.3883 




a 2 [GeV" 1 ] 




0.1755 


0.1575 


0.3246 


Si 


Ql 1 [GeV 2 ] 




5.2645 


4.5169 


3.9774 




s El [GeV 2 ] 




1.14° 


1.3038 


1.14° 




Gb 2 








-0.6520 




QO 








-0.8929 




a 2 [GeV" 1 ] 








1.7729 




Q% 2 [GeV 2 ] 








2.4634 




s E2 [GeV 2 ] 








1.14° 




so [GeV 2 ] 


1.14° 


1.14° 


1.14° 




Qo [GeV 2 ] 




0.4089 


0.4580 


0.9998 




Ql [GeV 2 ] 




3.1709 


2.5180 


1.8926 




N [GeV" 2 ] 


0.0408 


0.0655 


0.0567 




"Kd.o.f. 


12.92 


4.6886 


1.3005 



Table 1. Parameters of the fits. The symbol ° refers to the fixed parameters. 
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To start with, we made a fit by keeping some of the parameters fixed, close to their physical values, particularly 
those of the Regge trajectories and of the photoproduction amplitudes. Also, a single-term background was used. 
The resulting fit (fit 1) is shown in Table 1. Next (fit 2) some of the parameters of the Regge trajectory were varied. 
Consequently the x 2 w & s improved, although still remaining unsatisfactory. Finally, we let all the parameters vary (fit 
3) with the result reported in the Table. Fit 3 is good, with Xd.o.f. — 1-30. It is worth mentioning that a comparison 
with a similar fit performed in [2£j leading to Xd.o.f. — 9.4 needs care, since in [23] only one term in the background 
was included, the helicity amplitudes were kept constant and the dataset used included both data from and [l7| . 

To show the progress in the fits, we plot against the experimental data the structure functions for four different 
values of Q 2 with the parameters from three different fits - see Fig. . 

Having fitted the parameters (from now on we will use parameters of fit 3) , we can now proceed to further calculations 
(moments) and analyses (duality relations) of the model. 

V. MOMENTS 

We have calculated the moments of the structure functions using the explicit expressions and parameters fitted in 
the previous section. These moments can be used, in particular, to estimate the role of the non-perturbative effects 
(higher twists). 




10 L — — ■ — 1 — 1 — 

0.01 0.1 1 10 100 

Q 2 



FIG. 4: Nachtmann moments, M n for n = 2, 4, 6, 8. The plot compares the moments calculated from the Regge-dual with 
those extracted from the data and reported in \V!\ (inelastic part). 

From the operator product expansion (for a comprehensive review see e.g. [2lJ) the moments M n {Q 2 ) of F 2 are 
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defined as 



2 v (J-2)/2 

Mn(Q 2 ) = i?nj(^, Q Z )O n ,J^) [ ^ ) , (23) 



J=2k 



Q' 2 



where k = 1, 2, /i is a factorization scale, O n j([i) is the reduced matrix element of the local operators with definite 
spin n and twist J, related to the non-perturbative structure of the target, E n j(/jL, Q 2 ) is a dimensionless coefficients 
related to the small distance behaviour. 

The leading twist term r = 2 is well established in pQCD, while higher twists are indicators of the non-perturbative 
and confining effects. In order to study the higher twists, it is essential to have a complete knowledge of the F 2 covering 
the entire x— range for each fixed Q 2 . Higher twists can be well established only with higher moments (n > 2), 
meanwhile for M 2 their contribution is small even at Q 2 ~ 1 GeV 2 . Therefore the most interesting kinematical region 
lies between and 5 GeV 2 and large values of x, where the higher moments dominate. The JLab data and relevant 
calculations in 0] cover most of this region. 

In the present section we evaluate the Nachtmann (N) and Cornwall-Norton (CN) moments within our Regge-dual 
model and compare them with the data of the CLAS colaboration 01 as wen as with those from ref . (3^ . 

The relevant moments arc defined as 

Mi{Q 2 )= I dxp^Fz^Q 2 ) (24) 

where 



n 



v I (x)-l ) fOT I=N 

Pn( j ~ V- 2 , forI = CN 



"3 + 3(n + l)r + n(n + 2)r 2 " 
. (n + 2)(n + 3) J ' 

V/1 + 4MV7Q 2 , 
2x/(l+r) . 

Please note that in our calculations the elastic part of the SF (for x = 1) was not taken into account (see section III.G 
in Ref. [13). 

It is a relatively simple task to obtain the moments by using the existing numerical integration methods. We have 
used the parameters of fit 3 from Table 1. In Fig. 0]we plot the Nachtmann moments for n = 2, 4, 6, 8 together 
with the results from ^ n Fig- EI the calculated N- and CN- moments are compared with those from . On this 
second set of figures the errors in the momenta are not displayed; according to |32| they should be less than 5%. 

As seen from the figures, the agreement between our model and the data is quite good in the region Q 2 < 5 GeV 2 , 
where the SFs were fitted to the data. The discrepancies increase with Q 2 , away from the measurements. 



P(x,Q 2 ) = 
r = 



VI. DUALITY RATIO 



In this section we check the validity of the parton-hadron duality for our Regge-dual model by calculating the 
so-called 'duality ratio' 

I(Q 2 ) = J^- (25) 

^Scaling 

where 



Iscalingi^Q ) — / ds Fr^ 

J Smin 

I R es(Q 2 ) = f maX d S F 2 ReS , 



and we have fixed the lower integration limit s min — s , varying the upper limit s max equal 5 GeV 2 and 10 GeV 2 . 
These limits imply "global duality", i.e. a relation averaged over some interval in s (contrary to the so-called "local 
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FIG. 5: Nachtmann moments, , and Cornwall-Norton moments, M„ N , for n = 2, 4, 6, 8. These plots show the 
comparison between the moments evaluated according to our Regge-dual model and the values of the moments extracted from 
the electron-proton scattering data reported in (inelastic part). 



duality", assumed to hold at each resonance position). For fixed Q 2 the integration variable can be either s (as 
in our case), x or any of its modifications (a/, £, ...) with properly scaled integration limits. The difference may be 
noticeable at small values of Q 2 due to the target mass corrections (for details see e.g. 0])- These effects are typically 
non-perturbative and, apart from the choice of the variables, depend on detail of the model. 



In choosing the smooth "scaling curve" F^ a 19 (actually, it contains scaling violation, in accord with the DGLAP 
evolution) we rely on a model developed in |33| and based on a soft non-perturbative Regge pole input with subsequent 
evolution in Q 2 , calculated [33| from the DGLAP equation. 
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The function F 2 fles is our SF with the parameters of fit 3 (see Table 1). The results of the calculations for different 
values of s max are shown in Fig. [|J] 




Q 2 , GeV 2 

FIG. 6: Global parton-hadron duality test for different values of s ma x 

Given the available variety and flexibility of the existing parameterizations for the SFs (see Sec. Ill) we do no 
attribute too much importance to the above duality test. Its validity or failure to large extent may be caused by 
accidental interplay of the details of different parameterizations. By this we do not intend to raise doubts about the 
very concept of parton-hadron duality. Moreover, in our opinion, explicit realization of this concept, similar to the 
Veneziano model, should exist and looked for. Work in this direction is in progress. 

VII. CONCLUSIONS 

The main objective of the present study is a phenomenological analyses of the CLAS data in a model within 
the analytical S matrix approach, complementary to approaches based on pQCD. This analyses, as well as similar 
attempts show that achieving good fits (with low x 2 ) to the data is a highly nontrivial task by itself. The origin 
of this difficulty is the large number and high statistics of the data and poor understanding of the non-perturbative 
dynamics, typical of the kinematical region where data are collected. 

As repeatedly stressed, our approach does not compete with QCD; it is aimed to be complementary to QCD in the 
non-perturbative domain. The main virtue of our Regge-dual approach is its generality: potentially, it can be used for 
any value of its kinematical variable. From this point of view, of special interest is the possibility to link low-energy, 
resonance physics (and the JLab data) with the high-energy (or low x) physics (from HERA) by "Veneziano duality" 
(apart from parton-hadron duality), inherent in the model. 

The price for such generality is the available freedom or flexibility of the model. It can be, however, further limited 
by comparison with other models, pQCD calculations and the data. In particular, 

1. Realistic parameterizations for baryonic trajectories, satisfying the theoretical constraints yet fitting the data, 
should be further elaborated. Work in this direction is in progress. 

2. The separation of resonances from back grou nd is model-dependent. Our parameterization of the background 
differs from that introduced long ago (see e.g. [lS^ and used in all subsequent papers (e.g. Its non-orthodox 
motivation comes from dual analytical models. At the same time, fits to the data produce (see Sec. IV) a negative 
sign in front of the second term of the background, similar to the "orthodox" models (e.g. 0>E1)- 

3. The present Regge-dual approach generalizes the concept of transition form factors, continuous in spin. Moreover, 
higher spin resonance excitations are accompanied by higher powers of the relevant transition form factor, and since 
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the Regge trajectories imply an analytic continuation in spin, the same applies for the transition form factors. 

On the whole, the revival of the analytical methods, namely the study of various Riemann sheets of the scattering 
amplitude in the resonance region (for a recent interesting approach along these lines see (3J|), and its combinations 
with the parton model and QCD is a promising new development in the strong interaction theory, that may shed new 
light on the confinement problem. 

In estimating the predictive power (or flexibility) of the present model, we notice that the number of the free 
parameters here (23) is comparable or smaller to that in similar fits. For example I. Niculescu [35| uses 30 fitting 
parameters. The virtue of the present Regge-dual approach is the possibility to extend the model using the same 
set of the parameters to the small x domain, treated in Refs. |l4t Vwi l36|. Matching the large-x (Jlab) and small-x 
(HERA) kinematical regions will remove or at least reduce substantially the number of the free parameters and the 
constrain the flexibility of the model. The realization of this ambitious goal, already discussed in Refs. [hH ITsl l3t| . 
will depend on the right choice of the Q 2 dependence or, alternatively, the correct off mass shell continuation of the 
dual amplitude. In the present paper Q 2 — dependence was introduced in the resonance region via the transition form 
factors. 

To conclude, let us once more emphasize that the Regge-dual approach to DIS and DVCS to large extent is 
complemental to the conventional one, based on the presence of a hard scale, when Q (or a mass M) is large and the 
amplitude is calculable up to corrections of 1/Q times logarithms of Q 2 . In this case hard-scattering factorization can 
be applied for any x, small or not. 

In the standard approach the generalization of DIS structure functions to the DVCS amplitude can be illustrated 
(37) by the following sequence of transitions 

F 2 ~ 1m A(j*p — * 7*p) — * Jm A(j*p — > yp)t=o — ► A(^*p — > jp)t=o — * A(~/*p — > jp) . (26) 

In phenomenological approaches, t— dependence usually is introduced by simply multiplying the forward scattering 
amplitude by arbitrary exponential e Bt , incompatible with the shrinkage of the cone. A consistent, non-factorizable 
form of the t— dependence was discussed and derived within pQCD in a recent interesting paper by Freund Q. 

In the Regge-dual approach, on the other hand, the above sequence can be inverted: on starts with a complex, 
t— dependent DVCS amplitude that can be reduced to the DIS structure function F% by taking its imaginary part, 
setting t = and equating the two photon momenta. This approach does not require the presence of any hard scale, 
such as large photon momenta. The external photons are assumed to couple to the proton by vector dominance (or 
generalized vector dominance |l3j|). In this sense this approach is typically "non-perturbative" . Partons (quarks and 
gluons) are not present explicitly but rather implicitly, manifest in the scaling behavior of the amplitude for large s, t 
and/or Q 2 , as well as in the values of the parameters (e.g. quark counting). The link between the scaling behavior of 
the analytic and quark models is a very interesting but still open problem. It was approached in a number of papers, 
e.g. in |25| . where the large angle scaling behavior in a dual model was achieved by using Regge trajectories with 
logarithmic asymptotic behavior. 

Although ours is a typically "soft" approach, the quark structure, small-distance effects, etc are also present their 
due to the use of nonlinear Regge trajectories. In particular, the asymptotic logarithmic behavior of these trajectories 
could mimic hard scattering, quark counting etc. [2^, [33. These effects are not factorized, as in the standard 
approach of Q] and in most of the related papers, but are continuous, i.e. the transition from "hard" (perturbative) 
to "soft" (non-perturbative) dynamics occurs smoothly, according to the properties of dual analytical models. The 
correspondence between the "hard" sector of this dual model and pQCD (or the quark model) (see e.g. 24, 38] ) is an 
interesting problem, meriting further studies. 
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APPENDIX A: POLE DECOMPOSITION OF THE DUAL AMPLITUDE AND THE BREIT-WIGNER 

FORMULA 



In the vicinity of a resonance, the nucleon structure function can be can be written in a factorized form [29j : 

F 2 (x, Q 2 ) = -^-MW 2 - M 2 ) x P F (Q 2 ) (Al) 
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where Pp(Q 2 ) stands for some power of the nucleon (transition) form factor: this power is two in the standard 
approach, as e.g. in ref. 0, 0, 0, |2(|, but varies (rises) with the resonances spin in the present Regge-dual 
approach; v = 2^2 = JaL-. (p j s the four-dimensional momentum of the nucleon, q is the four-dimensional momentum 
of photon, see Fig. and M is the mass of the resonance. 

This formula determines the contribution of a single, infinitely narrow resonance to nucleon structure functions. For 
a wide resonance, if we replace the delta-function 8(W 2 — M 2 ) in the above expression by the familiar Breit-Wigner 
formula 



1 Mr 



7T (W 2 — M 2 ) 2 + M 2 T 2 ' 
where V is the resonance width, what leads to the following expression (2(|: 



(A2) 



2 2m 2 x 1 Mr 2 

F 2 {X,Q )= 1+4m2x2/Q2 -( m 2 + Q2 (1/x _ 1) _ M 2 ) 2 +M 2 r 2 * P f{Q )■ (A3) 

Now let us compare this expression with our eq. 11' I'D : 

F ^ - taffUm^-^t^ x Ff{q2) • (A4) 

Expanding the Regge trajectory near a resonance: IZej w rij + {7?.e ay} (s — M 2 ) = rij + {IZe ay} (m 2 + Q 2 (l/x — 
1) — M 2 ) and introducing the notation: T = j^~yM ' we ^ e ex P ress i on: 

2 Q 2 (l — x) N MY 2 

F2 ^' Q ' = 47ra(l + 4m 2 x 2 /Q 2 ) {^e a,}' (m 2 + Q 2 (l/x - 1) - M 2 ) 2 + M 2 T 2 X Pf ^ Q ^ ' (A5) 

Notice that Q 2 (l — x) — (s — m 2 )a; w [M 2 — m 2 )x in the vicinity of the resonance and therefore eqs. (|A3|) and (|A5J) 
are approximately the same for 

N = gm^ge a }' 
(M^ — m^) 

The obtained value for the normalization coefficient is approximately (for M = \/2 m and {7?.e aj} = 1 GeV~ 2 ) 
N w 8a = 0.058 GeV -2 , in agreement with the results of the fit (see Table 1). 
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FIG. 7: Structure function F 2 (x) for Q 2 = 0.225 - 0.925 GeV 2 . Data are from |17|. whereas the straight line is the prediction 
of our dual model. 




FIG. 8: Structure function F 2 (x) for Q 2 = 0.975 - 1.675 GeV 2 . 




FIG. 9: Structure function F 2 (x) for Q 2 = 1.725 - 2.425 GeV 2 . 




FIG. 10: Structure function F 2 (x) for Q 2 = 2.475 - 3.175 GeV 2 . 




FIG. 11: Structure function F 2 (x) for Q 2 = 3.225 - 3.925 GeV 2 . 




FIG. 12: Structure function F 2 (x) for Q 2 = 3.975 - 4.675 GeV 2 . 



